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Abstract. - We consider a discrete quantum system coupled to a finite bath, which may 
^ ' consist of only one particle, in contrast to the standard baths which usually consist of continua 

I of oscillators, spins, etc. We find that such finite baths may nevertheless equilibrate the system 

. though not necessarily in the way predicted by standard open system techniques. This behavior 

' results regardless of the initial state being correlated or not. 

in 

■ Due to the linearity of the Schrodinger equation concepts like ergodicity or mixing are 
_ strictly speaking absent in quantum mechanics. Hence the tendency towards equilibrium is 

not easy to explain. However, except for some ideas [1,2] the approaches to thermalization 
\ in the quantum domain seem to be centered around the idea of a thermostat, i.e., some 
^ ■ environmental quantum system (bath, reservoir), enforcing equilibrium upon the considered 
^ \ system. Usually it is assumed that this bath's classical analogon contains an infinite number 
O^' of decoupled degrees of freedom. 

Theories addressing such scenarios are the projection operator techniques (time-convolu- 
tionless, Nakajima Zwanzig), the Born approximation (BA) [3] and the path integral technique 
. ^ (Feynman Vernon [4]). The projection operator techniques are exact if all orders of the system- 

■ bath interaction strength are taken into account which is practically unfeasible. However, 
assuming weak interactions and accordingly truncating at leading order in the interaction 
strength (BA) produces an exponential relaxation behavior (c.f. [5,6]) whenever the bath 
consists of an continuum of oscillators, spins, etc. The origin of statistical dynamics is routinely 
based on this scheme, if it breaks down no exponential thermalization can a priori be expected. 

We find that this scheme breaks down (i.e. the BA produces wrong results) if the bath 
features a special spectral structure which cannot arise from an uncoupled multitude of sub- 
systems or modes (see below). We refer to this type of bath as finite bath. This holds true 
even and especially in the limit of weak coupling and arbitrarily dense bath spectra. 

Nevertheless a statistical relaxation behavior can be induced by finite baths. It simply is 
not the behavior predicted by the BA. Thus the principles of statistical mechanics in some 
sense apply below the infinite particle number limit and beyond the BA. 

(*) jgenmier@uos.de 
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Fig. 1 - Two-level system coupled to a finite bath. The environment's spectrum deviates significantly 
from those of baths with infinitely many decoupled degrees of freedom 



This also supports the concept of systems being driven towards equilibrium through in- 
creasing correlations with their baths [7-11] rather than the idea of system and bath remaining 
factorizable, which is often attributed to the BA [3,4]. 

The model we analyze is characterized (for simplicity) by a two level system (S or "spin") 
coupled two a many level system (B) consisting of two relevant bands featuring the same 
width and equidistant level spacing (see Fig.^. So this may be viewed as a spin coupled to a 
single molecule, a one particle quantum dot, an atom or simply a single harmonic oscillator. 
Note that the spin, unlike in typical oscillator baths or the Jaynes-Cummings Model, is not 
in resonance with the environments level spacing but with the energy distance between the 
bands. There are two principal differences of such an finite environment level scheme from 
the level scheme of, say, a standard oscillator bath, i) The total amount of levels within a 
band may be finite, ii) Even more important, from, e.g., the ground state of a standard bath 
there are infinitely many resonant transitions to the "one-excitation-states" of the bath. But 
from all those, the "back- transitions" lead to only one ground state. Thus the relevant bands 
of any infinite bath would consist of only one state in the lower and infinitely many states in 
the upper band. And it will turn out to be that limit in which the standard methods produce 
correct results. 

A finite bath which cannot be decomposed in uncoupled subunits any further may, however, 
feature arbitrary numbers of states in both bands. Note that in our model there is no notion 
of the environment being chaotic in itself. Due to the considered type of interaction the full 
system might be termed chaotic, as will become clear below (for a treatment of finite baths 
under a different perspective, see [11,12]). The Hamiltonian of the model in the Schrodinger 
picture reads H = Hq + V, Hq representing the uncoupled system and V the interaction: 

Ho = A£;(T, + V-fni|ni)(ni| -f V(A£;-H n2)|n2)(n2| 

iVi ^ — ' iVo 

V = XJ2 C(ni,n2)<7+|7ii)(n2|+h.c. (1) 

"1 ,"2 

Here the Pauli matrices refer to S, 711(71,2) denotes the 77,'th energy eigenstate within the 
lower(upper) band of B and h.c. stands for the hermitian conjugate the previous sum. For 
the example at hand we chose AE = 25u, 6e — 0.5u, Ni — N2 — 500, A = 5 • 10^'^u, u being 
some arbitrary energy unit. The real and imaginary parts of the C"s are randomly (Gaussian) 
distributed numbers with mean zero and normalized to X)ni na "■2)l^/-^i-^2 = 1- This 

interaction type has been chosen in order to keep the model as general and free from pecu- 
liarities as possible. (In the fields of nuclear physics or quantum chaos random matrices are 
routinely used to model unknown interaction potentials. We do, however, analyze the dynam- 
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Fig. 2 - Evolution of the excitation probability (product initial state) (a) and the off-diagonal element 
(correlated initial state) (b) of a spin coupled to a finite bath. Regardless of the baths finiteness the 
evolution exhibits statistical thermalization respectively coherence vanishes, i.e., the system evolves 
towards a maximum equilibrium entropy. For the diagonal elements, although, different from what 
the Born-approximation predicts. 



ics generated by one single interaction, not the average dynamics of an Gaussian ensemble of 
interaction matrices.) 

We firstly analyze the decay behavior of two different pure product initial states: The bath- 
part of both initial states is a pure state that only occupies the lower band but is, apart from 
that, chosen at random. Apart from its pureness only with respect to occupation numbers, 
B's initial state can be considered an approximation to a Gibbs state with Se <C fcTe ^ AE 
(in the example at hand, e.g., kT^ « 5u). For small Se the temperature may be arbitrarily 
small. Initially, the system S is firstly chosen to be completely in its excited state (this initial 
state is indicated by the black dots in Fig.^l and, secondly, in a 50:50 superposition of ground 
and excited state. The probability (density matrix element pii(t)) to find the system excited 



as produced by the first initial state is shown in Fig. 2(a) Since the first initial state does 
not contain any off-diagonal elements, we find |poiP ~ for all times. This is different for the 
second initial state investigated in Fig. 2(b) it starts with IpoiP = 0.25 and is thus well suited 
to study the decay of the coherence. (The diagonal elements of the second state are already 
at their equilibrium value (pii(O) = 0.5) in the beginning and exhibit no further change.) 

By numerically solving the time-dependent Schrodinger equation for the full model's pure 
state |\E'(t)) we find for the reduced state of the system p{t) = TrB{|vE'(t))(vE'(t)|}, an exponen- 
tial decay, up to some fluctuations. (For the baths initial state being a real mixed Gibbs state 
one can even expect fluctuations to be smaller, since fluctuations corresponding to various 
pure addends of the Gibbs state will partially cancel each other.) Thus the Schrodingerian 
dynamics yield a local behavior that might be described as statistical, controlled by some 
transition rate. 

The full model is Markovian in the sense that bath correlations decay much faster than the 
system relaxes, concretely bath correlations decay on a time scale of Tc ^ h/6e = 2 (all times 



given in units of h/u), whereas the system relaxes on a timescale Tr « 640 (cf. Fig.'s 2(a) (b)). 
Since the whole system is finite there is a finite (quasi) recurrence time. But due to the 
incommensurability of the full system's frequencies it appears to be > 10®, i.e., orders of 
magnitudes larger than the relaxation time of S. (The special example at hand features, due 
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to the environments equidistant level spacing, a recurrence time for the bath correlations of 
approximately 6 • 10'^ but that does not induce a recurrence in S.) 

Although the model is Markovian in the above sense and its relaxation appears locally 
statistical, S's excitation probability deviates significantly from what the BA predicts (cf. Fig 



2(a) I: The beginning is described correctly, but rather than ending up at T = Tb as the 
BA predicts for thermal environment states [3], S ends up at T = oo, i.e., equal occupation 
probabilities for both levels. Furthermore a condition often attributed to the BA, namely that 
S and B remain unentangled, is not fulfilled: When S has reached equilibrium the full system 
is in a superposition of |S in the excited state (g) B in the lower band) and |S in the ground state 
(8) B in the upper band). This is a maximum entangled state with two orthogonal addends, 
one of which features a bath population corresponding to Tb ~ 0, the other a bath population 
inversion, i.e., even a negative bath temperature. These findings contradict the concept of 
factorizability but are in accord with a result from [13] claiming that an evolution towards 
local equilibrium is always accompanied by an increase of system-bath correlations. However, 
the off-diagonal element evolution coincides with the behavior predicted by the BA. Thus in 
spite of the systems finiteness and the reversibility of the underlying Schrodinger equation 



S evolves towards maximum local von Neumann entropy (see Fig. 2(b) i which supports the 
concepts of [8]. 

We now very roughly (and rather incompletely) outline the Hilbert space Average Method 
(HAM) which explains the behavior of this model-type. HAM is not limited to two level 
systems, the example has just been chosen for simplicity. (For a detailed description of HAM, 
see [14, 15]). We start by considering short time steps of the evolution of S's density matrix p. 
With a truncated (second order) Dyson series for \^{t + At)) k, D{t, At)\'if{t)), one formally 
gets 

p{t + At)^TTB{D\^it)){^{t)\D^. (2) 

If the right hand side of Q (r.h.s.) was only a function of S's local density matrix p, one could 
set up an iterative scheme for the local dynamics of S. But one finds that the r.h.s. depends 
on all the details of the full model's state |^'(i)) and explicitly on the absolute time t, thus no 
autonomous iterative sub-dynamics for S can be inferred directly. However, computing the 
Hilbert space average (|- • •]) of the r.h.s., i.e., the average over an adequate set of full system 
states |<i>) sharing (some) crucial quantities with \^{t)), i.e., ($|^|$) = (^'(i)|A|*(t)) with A 
being the full model's energy, the local energy of S or the coherence of S , yields under specific 
conditions on the model (see below) 

iTTB{m){^\Dm « m + At {cm +im) , m 

where £, T are linear super-operators. (Note that the majority of the states belonging to 
the above set is correlated, i.e., no factorization is implied.) Replacing the r.h.s. of ||2J) by its 
Hilbert space average would thus yield an autonomous iteration scheme. But is that justified? 
It is justified, whenever Q does not only hold for the Hilbert space average over all |<&)'s 
but also (approximately) for the majority of all |<&)'s belonging to the above set individually. 
If this is the case the evolution of S tends to be independent of the full models state and 
basically controlled by the local state of S. This tendency increases with B featuring increas- 



ing numbers of eigenstates as can be shown theoretically [15] or numerically (cf. Fig. 3(b) I. 
Nevertheless, this replacement does not represent an approximation with a computable error, 
it only represents a best unbiased guess for the evolution of S. This "best guess" structure of 
HAM accounts for the non-statistical character of the under-laying dynamics, i.e., allows for 
a non-statistical behavior of S in some (rare) cases. 
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Performing the above replacement and taking Ai to zero yields for the system-type at 
hand (but independent of the concrete interaction) the following master equation scheme: 

, , , , , 2tt\^N2 2tt\^Ni 

Pii(i) = -(-Rio + -Roi)pii(i + -RioPii(O) , i?oi = — , itio = — rr 

noe hoe 

Poi(i) = {tAE/h - i?oi/2)poi(t)- (4) 

The equilibrium value of S's excitation probability is given by pii(oo) = Ni/{Ni + N2). Thus 
only if N2 ^ A^i (infinite bath) the BA produces correct results. Otherwise BA must fail for 
the reduced sub-dynamics of S are not even Markovian now in the sense of the equilibrium 
state being independent of the initial state. Note, however that it is not the finite density of 
states that causes the break down of the BA, since the BA produces wrong results even for 
Ni,N2 ^ cxD as long as the above condition is not met. Since boldly calculating transition 
probabilities according to Fermi's Golden Rule, would have produced the same rates R the 
applicability of the above scheme implies the applicability of a random phase approximation 
in the sense of, say, Peierls [16]. 

Such an approximation cannot always hold since the Schrodinger equation is completely 
non-random. Thus, as necessary conditions on the model parameters for statistically appear- 
ing behavior of S we find (from theory and numerics) 

2a£ > 1, X'^ « 1, (5) 

where N refers to the band with the larger state density. Those conditions enforce Markovicity 
in the above sense, and exclude the infinitely weak coupling limit, A ^ as long as B's state 
density is finite. Neither the equidistant level spacing of B nor the Gaussian distribution of 
the interaction matrix elements are indispensable. A condition on the level structure of B is 
that the number of states within an interval of / « 6e/l0 does not depend much on where 
within the band the interval is chosen. A similar condition restricts the interaction: The sum 
of |Cp's corresponding to transitions from one state of a band to states within an interval / 
of the other band should not depend much on the position of the interval. 

Since HAM is just a "best guess theory" the exact evolution follows its predictions with 
different accuracies for different initial states, even if all conditions on the model are fulfilled. 
To analyze this for, say pii(t), we introduce D^, being the time-averaged quadratic deviation 
of HAM from the exact (Schrodinger) result 



l£ {pfnt)-pir\t)ydt. (6) 

Thus D is a measure of the deviations from a predicted behavior. The results of the investi- 



D" = 

T 



gation for our model (Fig.QJ are condensed in the histogram (Fig. 3(a) t = 2000). The set 
of respective initial states is characterized by a probability of 3/4 for |S in its excited state (8) 
B in its lower band) and 1/4 for |S in its ground state (g) B in its upper band). Within these 
restrictions the initial states are uniformly distributed in the corresponding Hilbert subspace. 
Since all of them are correlated the application of a product projection operator technique 



would practically be unfeasible. However, as Fig. 3(a) shows, the vast majority of them follows 
the HAM prediction quite closely, although there is a typical fluctuation of D = • 10~^ 
which is small compared to the features of the predicted behavior (which are on the order of 
one), due to the finite size of the environment (cf. also fluctuations in Fig. |2(a ~ 



In Fig. 3(b) the dependence of on the number of states of B is displayed for N 



10,..., 800 ( one evolution for each environment size). At = 500 like used in the above 
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Fig. 3 - a) Deviation of the exact evolution of the spins excitation probability from the HAM predic- 
tion for a set of entangled initial states, b) Deviation of the exact evolution of the spins excitation 
probability from the HAM prediction for increasing number of states in the environment. 



accuracy investigation we find the same typical fluctuation, whereas for smaller environments 
the typical deviation is much bigger. We find that the squared deviation scales as 1 /N with the 
environment size, thus making HAM a reasonably reliable guess for many-state environments. 

What about the claim that reduced dynamics need not to be completely positive, once S 
and B are correlated [4, 17, 18]? In principle this holds. But regardless of their being correlated 
just a small fraction of all states from Fig. 3(a) shows significant deviations from HAM. Thus 
a smooth evolution towards equilibrium (though not necessarily of the Lindblad-type [19]) can 
typically be expected for the reduced dynamics, regardless of the initial state being correlated 
or not. 

What about reversibility and statistically appearing dynamics? The Schrodinger equation 
is time reversible and indeed there is no apparent time asymmetry in Fig. 2| For t > the 
behavior of the above system (solid line) is well described by the above rate equation scheme 
featuring an attractive fix-point. But for t < Q it would have to be described by a scheme 
with an repulsive fix-point! Thus, any of the full model states at i < is a paradigm for 
an initial state that does not yield statistical decay behavior eventhough the model typically 
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Fig. 4 - Evolution of the spins excitation probability before and after full excitation for two different 
sets of model parameters. 
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generates it. The dashed hne shows the behavior of a model hke the above one, only with 
(5e w and (in order to keep the timescales comparable) A = 10"**^. This clearly violates the 
second criterion of |(SJ and the model is no longer Markovian in the above mentioned sense. 
And indeed even for t > there is no exponential decay although the model features equally 
many states as the above one. Nevertheless the same equilibrium state as in the above case is 
reached (cf. also [20]) regardless of whether the system is propagated forwards or backwards 
in time. 

In essence we have shown that statistical relaxation may emerge directly from the Schro- 
dinger equation. This requires the respective system being coupled in an adequate way to a 
suitable environment. This environment must feature many eigenstates. There is, however, 
no minimum particle number limit. Thus the thermodynamic limit appears to be essentially 
controlled by the number of environmental eigenstates involved in the dynamics rather than 
by the number of environmental particles. This relaxation behavior results even for correlated 
initial states, nevertheless, standard open system methods may fail to produce the correct 
result. 

* * * 
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